Study of Bs — > r]c{J/i/j)D decay with perturbative QCD approach 

Junfeng Sun, Zhengjun Xiong, Yueling Yang, and Gongru Lu 
College of Physics and Electronical Engineering, 
Henan Normal University, Xinxiang 453007, China 

Abstract 

The T^-exchange process Bs — )• ric{J/tp)D is studied with the perturbative QCD approach. 
Three kinds of wave functions for Bs meson and two forms of wave functions for charmonium are 
considered. It is estimated that branching ratios for Bg — )• rjcD, rjcD, J/ipD, J/ipD decays are the 
order of 10^'', 10^^, 10^^, 10^^, respectively, where the largest uncertainty is from wave functions. 
There is a possibility for measuring these decay in the near future. 
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The first evidence for Bs production in e'^e anniliilation at tlie T(55') resonance was 



found by the CLEO collaboration jlj]. Bel 



e has accumulated 121 fb~^ data at the T{5S) 
resonance, including 7.1x10® -Bs-Bs pairs It is estimated that some 5.9x10*^ Bg mesons 
in the dataset of 5 ab^^ at the T{5S) resonance in a New Snowmass Year (about 10'^ seconds 
of actual annual running time jsj) will be collected at the high luminosity e'^e~ asymmetric 
SuperKEKB [4']. More and more B^ decays will have subjected the Standard Model and 
new physics to a series of increasingly stringent tests, through observables such as branching 
ratios, CP-violationg asymmetries and kinematic distributions. 

In the naive spectator model, the general properties of the Bg meson parallel those of the 
Bu^d mesons. The close correspondences between Bs and B^.d mesons allow for sensitive tests 
of hadronic models Hadronic B decays are complicated on account of strong interaction 
effects, meanwhile, they will have provided a great opportunity to study perturbative and 
nonperturbative QCD. For nonleptonic two-body B decays, the low-energy effective Hamil- 
tonian ansatz and the factorization hypothesis are commonly used. In recent years, several 
attractive methods have been proposed and widely used to evaluate the hadronic matrix ele- 
ments (where the local operators in the effective Hamiltonian are sandwiched between initial 
and final hadron states considered) based on an expansion in the QCD coupling constant 
as/ 71 and in the power Aqcd/'^q (where Aqcd and ttiq are the QCD characteristic scale and 
the mass of heavy quark Q, respectively), such as the QCD factorization jsl, perturbative 
QCD method (pQCD) [g], soft and coUinear effective theory (SCET) [3], etc. 

Using the operator product expansion and renormalization group equation, the low energy 
effective Hamiltonian for the Bg — )■ r]c{J/'ip)D decay can be written as jsj: 



T-ieff = -^V*iyus[Ci{baCa)v-A{u,3S0)v-A + C2(baCi3)v~A{upSa)v~A] 

G - - 

a Ua)v-A{cpSp)v-A + C2{haUp)v-A{Ci3Sa)v-A } + H.c., (1) 

where Gp is the Fermi coupling constant. The Cabibbo-Kobayashi-Maskawa (CKM) matrix 
factors V*iyus (or V*^ycs) and the Wilson coefficients Ci^2 describe the strengths of the local 
four-quark operators in the effective Hamiltonian. a and /3 are SU{?)) color indices. {qq')v-A 
= ^7^(1 — 75)0''. The Wilson coefficients, which incorporate the physics contributions from 
high scales, have been calculated to the next-to-leading order in the perturbation theory and 



evolved to a characteristic scale with the renormalization group equation 



. The essential 



problem obstructing the calculation of nonleptonic decay amplitudes is how to evaluate the 



hadronic matrix elements of the local operators properly and accurately. 

Based on the principle of color transparency 10| and factorization scheme, the phe- 
nomenological treatment of the hadronic matrix elements for the VT-exchange processes Bg 
—7- ric{J/ip)D is the same as that for pure annihilation topologies. Although the annihila- 
tion amplitude is formally power suppressed by Aqcd/'"^6 with the QCDF power counting 
conventions [.5|, its contribution is indispensable for realistic S-meson decays jo]. The com- 
prehensive analysis of Bu,d PP, PV decays without taking into account the annihilation 
contributions is generally of poor quality (P and V denote the pseudoscalar and vector 
mesons, respectively). Study of Bg — )• ric{J/i')D decay will help to improve our understand- 
ing of the annihilation effects. 

Analogous to the analysis fo. harf exclusive sca.e.ing a^pLudes Q. the had..on,c .a- 
trix element is commonly expressed as a convolution of scattering kernels with the universal 
wave functions of the participating hadrons , where nonperturbative dynamics either 
cancel or is absorbed into hadron wave functions (WFs). However, the annihilation effects 
in the coUinear approximation exhibit endpoint singularities (ES) for charmless mesonic B 
decays, displaying inconsistency of the QCDF formula [5|]. 

To deal with ES in convolution integrals, many attempts will have been made. 

(1) A phenomenological parameterization of ES in annihilation contributions is originally 
proposed by QCDF itself 5|, which ont only introduces uncertainties in the QCDF's pre- 
diction of observables, especially for annihilation dominated processes HI] , but also provide 
no constraint on magnitudes of strong phases relevant to CP violation. 

(2) ES is removed by separating the physics at different momentum scales using the zero- 
bin subtraction to avoid double counting of soft degrees of freedom in SCET [l^, while the 
imaginary part of the amplitude is also dropped at the leading power in a,^Aocp/rnh- 

(3) The infrared finite gluon propagator and running coupling constant [14!], or/and 
Cutkosky rules 3] for the quark propagators, are used to serve as a natural cutoff , which 
has already been applied to B decays into two mesons 1618]. However, it is claimed [l8] 
that different predictions on branching ratios can be obtained with different solutions of 
the Schwinger-Dyson equations for gluon propagator and coupling constant due to different 
truncations and approximations. 

(4) By keeping the parton transverse momentum kT, and employing the Sudakov factors 
to smear the double logarithm in QCD radiative corrections and to suppress the endpoint 
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contribution of hadron wave functions in small transverse momentum region, ES in collinear 
approximation is eliminated with the pQCD approach, and the strong phases are pertur- 
batively calculated {gI- An example is the recently renewed study on the pure annihilation 
decays Bg — t- tt+tt^ and B^^ — )■ K^K^ with the pQCD approach [19I] which are in good 
agreement with the CDF and LHCb measurements. 

Despite disputes as to which one of above treatments is more effective than others, we 
will study the Bg — ?■ rjc{J/ip)D decays with the pQCD approach to give a rough estimate 
of their branching ratios. Based on factorization, the typical expression for the decay 
amplitudes with the pQCD approach can be expressed as 

JdkC{t)H{k,t)<^{k)e-^ (2) 

where C, H, $, and are Wilson coefficient, hard-scattering kernel, hadron WFs, and 
Sudakov factor, respectively. The typical scale t depends on topology and process. For 
convenience, the kinematics variables are described by light cone coordinate. The momenta 
of the valence quarks and hadrons in the rest frame of the Bg meson are defined as follows: 

Pi = ^(1,1,0^), (3) 

P2 = {vt,V2.0±), (4) 
P3 = {%,vt,^±), (5) 
ki = XiPi + (0,0, fcix), (6) 

= 

62 = — (%+,-%~,0x), (8) 
1712 

where the subscript i = 1, 2, 3 refers to Bg, r]c{J/ip), D meson, fcj, fcj^, are the momen- 
tum, transverse momentum and longitudinal momentum fraction of light quark of meson, 
respectively. €2 is the longitudinal polarization vector of J/ip meson. In the rest frame of 
Bg meson, Ei is the energy of particle i, and p is the common momentum of final state. 

The basic input element in Eq.([2]) — WFs — is defined by the nonlocal bilinear quark 
operator matrix element 21 1. 



{Q\hMs,{z)\Bg{p,)) = -^ld'k,e-^'^-{[^^^^ (9) 



Pa 



1 



{J/^{p2)\-cMcp{zm = /d^fc2e+^^-Y2 mj/^<p^^{k2)+M'A''2) (10) 



{Tlc{p2)\cMcp{zm 



—I 



d^A;2 e 



+ik2-z 



{Dips)\c^{0)u^iz)\0) 



V2K 



{75 [mi ih) + m^,^(f)l (fca)] }^^, 



(11) 
(12) 



where A^c = 3 is the color number. n_ and n+ are null vectors, and rij^-ri- 
Here, the distribution amplitude of D meson given in |20| is used, 

fn 



1. 



1 + Cd{1-2x) 



where x 



1 — X. fr, is the decay constant. Cr, is a shape parameter. 



(13) 



For WFs of Tjc and J/ip mesons, and (p^ are twist-2; and 0^ are twist-3. They 



can be extracted from the Schrodinger state with dynamical potentials [21|, |22|. We will 
consider two kinds of WFs corresponding to harmonic-oscillator and Coulomb potentials. 
Their expressions are listed in j22|. One is the harmonic-oscillator (O) type 



05; (x, 6) 

0^(x,6) 



^iV,Wp{->[(^)%.^.^]}, 



A^4(x -x)^exp{ 



fvc T\TV - ( 



-XX 



TTlr 



-XX 



The other is the Coulomb (C) type 



iV;exp{ 



nir 



-XX 



2xx 

V 2xx ^ J J 
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— 4xx(l — f^) 


(x 


— xYxxmJj 




— 4xx(l — f^) 




{xxYmcb 




— 4xx(l — v'^) 




xxirich 


Jl 


— 4xx(l — v"^) 



Ki {rricb^Jl — 4xx(l — v^)) , 



Ki{mchJl — 4xx(l — f^)). 



Ki {iTLcb^Jl - 4xx(l - f2)) . 



(14) 
(15) 
(16) 
(17) 

(18) 
(19) 
(20) 
(21) 



where fjj^ and /^^ are decay constants, rric is the mass of c quark, h is the conjugate variable 
of the transverse momentum, u ~ 0.6 GeV and v ~ 0.3 [22] are shape parameters. iV^'* 



and N'"'^ are the normalization constants. The normalization conditions are 



dx 05;'*(x,O) 
dx 0;;^(x,O) 
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/j/V 



2^/2Nl' 

fvc 
2V2K' 



(22) 
(23) 



For WFs of Bs meson, there are two scalar compositions 0^ and 0^. Neglecting three 



particle amplitudes, the equation of motion for 0| is 23l. \24\ 



x) + X( 



0. 



(24) 



The relation of Eg. (1241) is sometimes referred to as "Wandzura-Wilczek relation" 25|. It is 
helpful in constraining models for WFs, which leads to <f)^{x) vanished at the endpoint and 
'~(x) = 0{1) for X — )■ 26|. Here we will investigate t 



iree models of WFs for Bs meson. 



The first one is the exponential (GN) type suggested in 



27|, i.e.. 



X, h) 
X, h) 



Bs 



IBs 



:A^(tj>^xexp 
^^GNexp 



xniBs 



XTUb 



1 



The second one is the Gaussian (KLS) type proposed in 



28 



(b wgn)^ 



29|, i.e.. 



,KLS- 



{x,b) 
(x,6) 



JBs j,t+ 2-2 

^^ATklsX X exp 



l/a;mB\2 1 



2 ^ wkls 



u 



Bs 



:iVKLsexp 



UJ- 



KLS^ 



exp 



KLS*- 



1 /xmBs\'^ 

2 ^ wkls 



(25) 
(26) 

(27) 



+ v^ms.c^KLsErf (-=-^) + Ckls}, (28) 
where the constant Ckls is chosen so that 0§^^~(1, 6) = 0. The third one is the KKQT type 



derived by QCD equation of motion and heavy-quark symmetry constraint 

fBs 2x 



30| 



I.e. 



iKKQT- 
Bs 



X, h) 
X, h) 



2^2W^J^ 



fi 



KKQT 

22/ 



2y2iVcW;^ 



KKQT 



e{y)Jo{mBsby/xy 
e{y)Jo{mBshy/xy 



(29) 
(30) 



where y = wrkqt - x. 

In Eq. fl25l — [30]) . Jb^ is the decay constant. Ui is the shape parameters. The normalization 
conditions is 

fss ^3^^ 



(/)r (x,0)dx — — 

Within the pQCD framework, the Feynman diagrams for Bs — t- rjcD decay are shown 
in FIGlH where (a) and (b) are non-factorizable topologies, (c) and (d) are factorizable 
topologies. After a straightforward calculation with the master formula of Eq.(l2]), the decay 
amplitudes can be written as follows 



A{Bs^VcD) = '^V:,V,s E A, 

V ^ i=a,b,c,d 



(32) 



The expressions of Ai are collected in APPENDIX. From the expressions, we can clearly see 
that both and contribute to the decay amplitudes. The branching ratios in the Bs 
meson rest frame can be written as: 

Bn{Bs^ri,D) = ^^\AiB,^VcD)\'', (33) 

where p is the center-of-mass momentum. 

The input parameters in our numerical calculation are collected in TABLE. [H If not 
specified explicitly, we shall take their central values as default input. 



Our study show that (1) contributions of 



which is consistent with pQCD's statement jo]. The interference between factorizable dia- 
grams FIGil] (c) and (d) is destructive, which is, by and large, in agreement with previous 



IGJH (a-c) can provide large strong phases. 



pQCD's estimate (for example, see |3j]). The strong phase difference between FIG|T] (c) 
and (d) is independent of model of WFs for Bg meson. The main contribution is from non- 
factorizable diagram FIGlT](b). (2) The dominant contribution is from Os/vr < 0.2 region, 
implying that despite the small phase space, these processes are calculated with peturbative 
theory due to hard gluon exchange, where the gluon virtuality scales as /c^ > (2mc)^. (3) 
There is very strong interference between contributions of WFs 0^ and 0^, between con- 
tributions of twist-2 and twist-3 WFs for rjciJ/ip) mesons. Contribution with only twist-3 
WFs for ric{J/ip) mesons (where twist-2 WFs is zero and twist-3 WFs is nonzero) is less 
than 30%. 

Our numerical results are shown in TABLE. |Tll where the first uncertainty comes from 
the WF shape parameter of Bg meson, i.e., cjgn = 0.45±0.10 GeV in Eg. fl25l— [26]) . wkls 
= 0.45±0.10 GeV in Eq.(l23H28]) and wkkqt = 0.25±0.10 in Eq. (I29H30]) ; the second un- 
certainty comes from the WF shape parameter of J/iplrjc) meson, i.e., u = 0.6±0.1 GeV 
in Eq. (jMl — ITTl) and v = 0.3±0.1 in Eq. (JT8l— I2T1) : the third uncertainty comes from the WF 
shape parameter of D meson, i.e., Cd = 0.7±0.1 in Eq.( !T3|) : the last uncertainty comes from 
the choice of hard scales (l±0.1)tj in Eq. flEllI — IE12|) . In addition, decay constants fa, fj/tp, 
fr)c-i fss bring some 10% uncertainty to the branching ratio. 

From the numbers in TABLE. [Tll we can clearly see (1) branching ratios are sensitive 
to the choice of shape parameter and model of hadronic WFs for Bg and rjclJ/ip) mesons, 
relative to the choice of hard scale. It is also found that all branching ratios decrease with the 
increasing shape parameter of hadronic WFs for Bg and ric{J/ip) mesons. (2) Due to CKM 
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factors |Vc*t,K^<l > lK*6^cs|, there is hierarchic structure B7l{B J / ipD) > BTZ^Bs^J/il^D) 
and B7l{Bs^T]cD) > BTZ{Bs^r]cD). Besides, uncertainty (~ 30%) from V*^Vcs is much 
larger than that (~ 5%) from V^^Kis- (3) Due to mj/^ > m.^^ and the orbital angular 
momentum -Z^j/^75(^) > L^jj^^j^^, the phase space for Bg — t- J/ipD{D) decay is tighter than 
that for Bs — )■ ricD{D) decay. With the same input, there are relations BTZ{Bs-^ricD) > 
Bn{Bs^J/tpD) and Bn{Bs-^r]cD) > Bn{Bs^J/i)D). (4) Branching ratios for fi, rjJD, 
rjcD, J/ipD, J/ipD decays are the order of 10"'', lO"*^, 10"*^, 10~^, respectively. 

The corresponding [/-spin process B^ — )■ ric{J/ip)D{D) has been studied |35l-l38j . In |35| . 
it is argued that if the intrinsic charm inside B meson is not much less than 1%, branching 
ratio for B^ipdcc) — > ric{J/'4')D decay is ~ 10"^, which is larger than the present experimental 
upper limit < 1.3x10"^ )9L Based on collinear factorization scenario, the B^ — ?■ ric.{J/ip)D 
decay is investigated in 36| with the approach for exclusive processes [12!], where narrow 



5-function like WFs are used. The small overlapping among WFs results in branching ratio 
being about 10"^ ~ 10"^ [36]. This issue is renewed in [37] with pQCD approach in the 
framework of factorization. By keeping the parton transverse momentum and taking the 
WFs for cc final states given in 2l|], it is found that branching ratio for Bd — )■ i]c{J/'ip)D 
decay is bout 10^^ ~ 10"'' 371]. Considering the final state interactions, branching ratio 
for Bd —7- J/tpD is estimated to be 10~^ ~ 10~^ [38]. The results in {sgI, [s^] have similar 
hierarchic structure due to kinematics and dynamics, i.e., BlZ^Bd^rjcD) > BlZ{Bd^J /ipD). 



The method used in our study is the same as 



371 ]. and similar WFs for ric{J/ip) is employed 



(in our study, the small relativistic corrections to the WFs are neglected and two types of 



37] . using the 



WFs are considered). A consistent estimation is obtained between ours and 

^^^^ Bn(B,^ccD) ^ KIVW ^ ^ ^ ^y^^ )■ 

It is well known that the pure annihilation process B^ — )■ tt+tt^ with branching ratio 



(9(10^'') [39] and pure leptonic rare decay Bg — )■ fi^fj^^ with branching ratio ~ O{10~^) 40] 
have recently been measured at hadron collider, due to the fact that there have accumu- 
lated much data and that detectors sitting at LHC and Tevatron colliders have excellent 
performance on the final charged particles. We believe that Bg — )■ r]c{J/ip)D decay could be 
accessible experimentally in the near future, because (1) their branching ratio is the same 
order as (sometimes larger than) that for Bg — t- tt+tt^, decays. (2) The final D meson 

can be tagged by charged kaon and/or pion, while tracks of both and vr^ are be clearly 
seen by sensitive detectors. Besides, signal of rjc{J/ip) meson is easily identified by its nar- 
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row peak in the invariant mass distribution. For example, the LHCb has measured many 
5, decay into final states containing a charmonium, such as Bg — ?■ J /ipK^K' [41[, J /ipK 
42|, J/^/o(980) y, J/'0/sr° M .... (3) More and more 5, data will be accumulated with 
the running of LHC and advancing SuperKEKB. It seems to exist a realistic possibility to 
study rare decays with branching ratio ~ (9(10^^). 

In summary, we study the pure weak annihilation process — )■ ric{J/ip)D decay with 
pQCD approach. ES disappear as expected by keeping the parton transverse momentum. 
The largest uncertainty in our result is mainly from QCD's dynamical property of hadron. 
Branching ratio for Bs — r]c{J/'ip)D decay depends strongly on model of WFs for Bg and 
rjclJ/ip) meson. There are some other uncertainties considered here, such as the high order 
corrections, the effects of final states interaction, and so on. Our estimate show that branch- 
ing ratios for Bg — ?■ t]cD, rj^D, J/tpD, J/ipD decays are the order of 10~^, 10^^, 10"^, 10^^, 
respectively. They could be measured in the near future. 
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Appendix A: The amplitudes for Bg — J/i]jD decay 



A{Bg^J/i^D) = '^V:,V^g Y: a 

i=a,b,c4 



iAa = — mi / dxi / dx2 / dxs / bidbi / b2db2 / db^ 
y2N Jo Jo Jo Jo Jo Jo 

X agita)Ciita)e-^^-'^^-^^Habia, Pa, bi, b2)Mx3, h)5ib2 - 63) 

mhT]^ + V2i]{xi - X2) + V2ml{xi - X3 



X { 0;^(xi,6i)0j;(x2,&2)% 

-(t)]^{xi,bi)(j)^{x2,b2)r]2 
+4>b{xi, 6i)0^(x2, b2)m2ms 

-05(a;i, 6i)0^(x2, b2)m2m3 



(Al) 



nihr]^ + V2r]{xi - X2) + V2ml{xi - X3) 

1 1 + 1+1 

TTib + -mixi - -^2^2 - 2:3 

m^^ -mxxx- —|^^^X2- -^3X3\], (A2) 
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327r^C^ rl rl roo rOQ roo 

lAh = — j^^nii / dxi / dx2 / dxs / 6id6i / 62d62 / dfes 
v2A^ -^0 -'o -'0 JO JO JO 



X as{tb)Ci{tb)e~ 



-S,b—S 



°Hab{a, /3b, bi, b2)(pD{x3, b3)S{b2 - 63) 



1 



1 



X {V2mip(l)^{x2,b2)[r]^(l)'^(xi,bi) +r]^ ^^(xi, 61)] (xi - x^) 

+m2m34>%{xi, bi)(plj,{x2, 62) 

-m2m30B(xi, 6i)0^(x2, 62) 



-mixi - -^7/2 a;2 - ;^^3"^3 



1 



1 



^rnixi - -j=ri2 X2 



1 

V2 



«A = j;^ — rnipfs. 



1 rl poo poo 

dx2 / dxs / 62d62 / b3dbs 
JO Jo Jo 



X a,(g[Ci(te) + ArC2(te)]e-^'^"^^i^cd(«,/9c,&2,&3)0D(a:3,63) 
X {(j)^{x2, 62) - - ml)x2 - 2m2m3X20^(x2, 62)}, 



iAd 



/I /"CO roo 

dx2 / dx3 / 62d62 / 
JO Jo Jo 



^—S^—Su 



b3db3 



X as{td)[C^{U) + NC2{U) 
X 0£,(a:3, 63)0j;(x2, b2)\ml + m3mc - (m^ - m2)x3}, 



Appendix B: The amplitudes for Bg J/ijjD decay 



i=a,b,c,d 



lA^ = ^^^==^mi / dxi I dx2 [ dx3 [ bidbi [ b2db2 [ db3 
y2N Jo Jo Jo Jo Jo Jo 

X as{ta)Ci{ta)e~^''-^^-^''Hab{a, 61, b2)M^3, h)S{b2 - bs) 
X { (I)'^{xi,bi)(j)^{x2,b2)rit V2mip{xi - X3) - rub^ 
bi)4>^{x2, 62)% V2mip{xi - X3) + mbr}2 



-(f)%{xi, bi)(f)lp{x2, b2)m2m3 



V3X3\ 



111 

+(t>B{xlMW^p{x2M)rn2rn3[-mxXl - -^2^1 - TTf'^a'^s]}, 



iA\y = ^^^==^mi / dxx I dx2 [ dx3 [ bidbi [ b2db2 [ dbs 
\/2N Jo Jo Jo Jo Jo Jo 
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X as{tb)Ci{tb)e 



-Sb—S^—Su 



Hab{oi, Pb, bi, h2)(t>D{xz, h)5{h2 - h) 



X {a/20^(,T2, 1)2) r7^0^(a;i, hi) - r]^ (Pb{xu h) vi^i - ^2) + ml{xi - X3) 
-m2m34>%{xi, 6i)0^(x2, ^2) 

+m2m3(f)B{xi, bi)(l)i{x2, 62) 



-mixi - -^i X2 - 



-mixi - X2 - X3 



Ac = -AdiBs^J/i^D), 
Ad = -Ac{Bs-^J/ijD) 

Appendix C: The amplitudes for Bg rjcD decay 



V i=a,b,c,d 



Aa, = ^^^==^mi / dxi [ dx2 [ dx^ [ bidbi [ 62(162 / dfea 
\/2N Jo Jo Jo Jo Jo Jo 



/2N 
X as{ta)Ci{ta)e 



—Sb—Stic—Sd 



Hab{a, Pa, bl, b2)(l)D{x3, &3)(^(&2 " ^s) 



X { 0^,(3^2, b2)(j)t{xi, hi)r]^ mb% + v^rjixi - X2) + V2ml{xi - X3) 
+(f)^^^{x2, b2)(f>B{xi, bi)ri2 [rribri^ + V2ri{xi - X2) + V2ml{xi - X3) 
+(t}'n,{^2,b2)(/)B{xi,bi)m2m3[mb + ^mixi - -^vt^^ - -^^t^s 

+0^,(2^2, b2)(f)B{xi, bi)m2m3 



1 1 _ 1 _ 

rUb + -miXi - -j=r]^ X2 - -^Vs X3 



/I pi pi POO POO POO 

dxi / dx2 / da;3 / bidbi / 62(162 / dfea 
Jo Jo Jo Jo Jo 



X a,(4)Ci(t,)e-^«-^''=-^«i7,5(a, 61, 62)0D(a;3, &3)5(&2 - ^3) 



X 



{V^(I>IS^2, ^2) [vt(l)B{xi, bl) + % (pBi^u bl)] [v{xi - X3) + ml{xi - X2) 
-0^^(x2, b2)(l)^{xi, bi)m2ms 



-miXi - -j=rij X2 - -^V3^3 



-(f>^r,ci^2, 62)0B(a;i, bi)m2m3 



-miXi - -j=r)2 X2 - -^V3 ^3 



Ac — 



N 



/I p1 poo POO 

dx2 / dx3 / 62d62 / bsdbs 
Jo Jo Jo 
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X { l^/^vt ^2 - mi\ [rit Vz<fnX^2, ^2) + m2rn3(/)^^(x2, 62)] 
+ [\/2772"x2 - mi] [jy^jys'^^el^s, ^2) + m2m30^^(x2, 62)] }, 



iV 



/I r\ roc roo 

Jo Jo Jo 



X asiU) [Ci{td) + NC2{td)]e-^^'^-^''Hcd{a, Pd, &3, b2)M^3, bs) 
X { [mi - V2r]^ X3] [r^^ r?^0^^(x2, 62) + m2m30^^(x2, 62) 
+ [mi - V2%X3 r]^r]^ 0^,(a;2, ^2) + m2m30^Ja;2, 62) 
-m3mcE2(l)l^{x2, 62) - 2m2mc-E30^^(x2, 62)}, 



Appendix D: The amplitudes for 5^ — > r/ci^ decay 



327r2CV 



V ^ i=a,b,c,d 



/I /•! POO roo roo 

dxi / da;2 / dx^ / 6id6i / 62d62 / d^s 
Jo Jo Jo Jo Jo 



X a,(t„)Ci(t„)e-^«-^''=-^^i/„,,(«, /3a, &2)0D(a;3, &3)'5(&2 - ^3) 
X { 0^^(x2, 62)05(2^1, &i)?7^ m^Ty^ + V2r](xi - X3) + V2ml(xi - X2) 
+05^^ (x2, b2)(i)B{xi, 61)773" [mbrj^ + \Pir]{xx - X3) + \/2m\{xx - X2) 

r 1 1 1-1 

+0^,(2:^2, b2)4>'^{xx, 6i)m2m3[m6 + -mi^i - -j^t ^2 - -^t^^ 
+0^^(x2, 62)0^(^1, bi)m2m3 



nib + -mixi 



= ^^^=^"^1 / dxi / da;2 / dxg / 6id6i / 62d62 / d63 
y2N Jo Jo Jo Jo Jo Jo 

X Q;,(ift)Ci(i6)e-^^-^''^-^^i/„ft(a, /J^, 61, 62)0D(a;3, 63)^(62 - 63) 

X {V20;j^(x2, 62) [nt (f^ti^u bi) + r7^0B(^i> ^1)] [^7(^1 - ^2) + m3(xi - X3) 



-(pt,,{x2, 62)05(0:1, 6i)m2m3 
-0^^(x2, b2)(l)B{xi, &i)m2m3 



-mixi - -^r]^ X2 - :^mX3 

1 _ 1 _ 1 _ 
-mia;i - -j=r]2 X2 - X3 
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Ac = -Ad{Bs^r]cD), 
Ad = -AciBs^VcD) 



(D4) 
(D5) 



Appendix E: Some formula in the decay amplitudes 

The 5*^ (5'^c,J/V' ^d) in the factor e^^^ (e~'^'"=''^/* and e~^°) is defined as 

Ssit) = s{x,pt,h) + 2f ^7,(/i), (El) 

J l/bi djJ, 

Sn.,J/i>{t) = s{x2pt, 62) + 2 /|/^^|^7,(/^), (E2) 
Suit) = s[x,p-,M) + 2f ^7.(/^), (E3) 

J 1/63 Ct/i 

where the quark anomalous dimension 7^ = — as/vr. The expression of s{Q^ h) is given in j^. 

H,,{a,PMM) = Ko{h^){e{h - h)HS,'\hV^)Jo{hV^) + {h^b2)}, (E4) 
H,,{a,P,b„bj) = H^'\b,V^){e{b,-b,)Ko{b,^)Io{b,^) + {b,^bj)}, (E5) 

The virtuahties of internal gluons (a) and quarks and the typical scale tj are defined 
as (where the subscript i = a, b, c, d corresponds to the FigJT]) 

a = xlml + xlml + 2x2X317, (E6) 
-(3a = ml{xi - X2)(xi - X3) - ml 

+ ml{x2 - Xi){x2 - X3) 

+ m3(x3 - xi)(x3 - X2), (E7) 
-(3b = mi(xi - X2)(xi - X3) 

+ m2(x2 - Xi)(x2 - X3) 

+ m3(x3 - xi)(x3 - X2), (E8) 

—(3c = m\ + x\m\ — X2{m\ + m\ — ml) , (E9) 

o 2 I 2 2 2 

—(3d = m]^+ x^mr^ — m^ 

- X3{ml-ml + ml), (ElO) 



ta(b) = max(A/a, yjA^, I/&1, 1/&2, 1/&3), (Ell) 
tc{d) = max ( Va, -^[5^, 1/62, 1/&3)- (E12) 
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TABLE I: input parameters for Bg — t- r]c{J/ip)D decay 



parameter 






value 


reference 


mass of Bs meson 




rriB, -- 


= 5366.77±0.24 MeV 


[9] 


mass of J/^l^ meson 




mj/^ = 


= 3096.916±0.011 MeV 


[9] 


mabo ui ineoUii 






= 2981. 0± 1.1 MeV 


fql 
I^J 


mass of D meson 




rriD = 


= 1864.86±0.13 MeV 


[9] 


mass of b quark 




mb 


= 4.18±0.03 GeV 


[9] 


mass of c quark 




nic 


= 1.275±0.025 GeV 


[9] 


lifetime of Bg meson 






= 1.497±0.015 ps 


[9] 


decay constant of D meson 




Id 


= 206.7±8.9 MeV 


[9] 


decay constant of J/^l^ meson 




fj/^ = 405±6 MeV 


m 


decay constant of r]c meson 




frjc 


= 394.7±2.4 MeV 


[32] 


decay constant of Bg meson 






= 227.6±5.0 MeV 


[33j 


Wolfenstein parameters 






A = 0.81ll°:°?i 








A = 


= 0.22535±0.00065 


[9] 








^ _ 101 +0.026 


[9] 








^ _ n O4r;+0.013 


[9] 
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TABLE II: Branching ratio for Bg — )• rjc{J/4^)D decay with different models of WFs. 



1.5< 
2.831J; 



GN 

94+0.61+0. 
61-0.34-0. 

02+3.56+0. 
74-1.73-0. 



13+0.05 
12-0.05 

09+0.25 
06-0.33 



KLS 



q no+0. 84+0. 96+0. 11+0.02 
"^•"-'^-0.91-0.62-0.08-0.03 

A 07+0.72+4.76+0.19+0.31 
' -0.89-2.41-0.16-0.42 



KKQT 



9 Q9+1. 63+0. 96+0. 08+0. 05 
^•^^-1.91-0.61-0.07-0.07 

A nn+1-58+4-88+0. 16+0.30 
^•"-'"-'-2.52-2.39-0.14-0.42 



Bn{Bs^J/'iljD)xW' 



4.34 



+1. 



7.08 



+0. 



37+0.27+0. 
40-0.16-0. 

65+1.45+0. 
55-0.98-0. 



08+0.31 
07-0.33 

07+0.36 
05-0.33 



6.28 



9.08 



+0.63+0.32+0.17+0.38 
-0.68-0.22-0.14-0.41 

+0.61+2.08+0.27+0.35 
-1.54-1.30-0.21-0.30 



7 nQ+O-32-l-l. 00+0. 12+0. 42 
' •"-'^-1. 84-0. 66-0. 09-0. 46 

1 c; 07+0-78-1-5.72+0.10+0.50 
10. U I --7,i5_3.48-0. 03-0.45 



4.09 
5.36 



+0. 
-0. 

-1. 



87+0.67+0. 
84-0.49-0. 

14+3.57+0. 
05-1.98-0. 



31+0.07 
27-0.09 

17+0.19 
11-0.24 



5.30 



7.22 



+0.55+0.89+0.36+0.09 
-0.68-0.64-0.31-0.11 

+0.63+5.02+0.16+0.24 
-0.98-2.74-0.07-0.31 



5.36 



+0.87+0.89+0.37+0.07 



-1.58-0.6 



-0.32-0.09 



7 9/1+1.14+4.77+0.19+0.21 
' ^^^-2. 74-2. 65-0. 10-0. 27 



Bn{Bs^r]cD)x 10^ 



9.48^ 



7.96 



+0. 



.08+0.21+0. 
.43-0.26-0. 

29+0.74+0. 
0.66-0. 



-0.86 



12+0.08 
10-0.00 

24+0.09 
22-0.01 



11.43 



-0.06+0.21+0.27+0.07 
-0.42-0.26-0.24-0.00 



12.46 



+0.02+0.58+0.22+0.09 
-1.95-0.55-0.19-0.00 



9.39 



+0.34+1.30+0.42+0.05 
-1.02-0.90-0.38-0.00 



11 nq+0. 15+2.31+0.33+0. 10 
li.yo_3 43_]^ g5_Q 28_o.oo 




FIG. 1: Feynman diagrams for Bg — ?• rjcD decay within the pQCD framework 
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